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VO ■ Abstract 
On 

t— i , 

We review briefly the fundamental equations of a semi-microscopic core- 
particle coupling method that makes no reference to an intrinsic system of 
<^ coordinates. We then demonstrate how an intrinsic system can be introduced 

qq ■ in the strong coupling limit so as to yield a completely equivalent formulation. 

It is emphasized that the conventional core-particle coupling calculation in- 
troduces a further approximation that avoids what has hitherto been the most 
time-consuming feature of the full theory, and that this approximation can 
be introduced either in the intrinsic system, the usual case, or in the labora- 
tory system, our preference. A new algorithm is described for the full theory 
that largely removes the difference in complexity between the two types of 
. calculation. Comparison of the full and approximate theories for some rep- 

! resentative cases provides a basis for the assessment of the accuracy of the 

traditional approach. We find that for well-deformed nuclei, e.g. 157 Gd and 
q ' 157 Tb, the core-coupling method and the full theory give similar results. 



O 



PACS number(s): 21.60.-n, 21.60.Ev, 21.10.-k, 21.10.Re 

*U ■ I- INTRODUCTION 

We have recently undertaken the task of revitalizing and extending a semi-microscopic 
theory of collective motion for odd nuclei that we shall refer to as the Kerman-Klein-Donau- 
Frauendorf (KKDF) model []I]-ffl. This model, aside from the elements discussed for the 
first time in the present paper, was introduced in close to its present form by Donau and 
Frauendorf whose work was in turn stimulated by an application [TTJ of the theory of 
collective motion developed by Kerman and Klein JTIj-Q. 

In the presentation of our work at seminars and conferences, one question that has 
invariably arisen is the connection between the KKDF model and the conventional core- 
particle coupling model, especially for deformed nuclei, to which our published applications 
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have so far been confined. Even if we widen the inquiry to the connection between the shell 
model and the core-particle model, we find that the literature on this subject is sparse. We 
are aware of only two publications that have been addressed specifically to this topic. The 
earlier of these papers |1| showed how all, then extant, core-particle coupling models could 
be understood as approximations to the work of Kerman and Klein. This paper appears 
to have gone completely unnoticed, since it is not quoted in the later work [16], which is 
devoted to the derivation of the strong coupling core-particle model from a schematic shell 
model. In the book by Ring and Schuck flT7| , which appeared betweentimes, the success 
of the strong coupling model in its domain of application is heralded but at the same time 
proclaimed a mystery. 

The main purposes of the present work are threefold. The first is to transform the 
Kerman-Klein equations from the "laboratory" system in which they are derived and con- 
veniently applied to the "intrinsic" system, when it makes sense to define such a system, 
as is done in the strong coupling core-particle model. The resulting theory is completely 
equivalent to the starting one and does not yet constitute the standard phenomenological 
model. A second purpose is to describe and implement the approximation that leads to 
the standard model. We describe in most detail how this may be done in the intrinsic sys- 
tem, the usual choice, but emphasize that the approximation may equally be defined in the 
laboratory system and that the latter approach has some advantages. 

The essential point here may be described as follows. In the physical situation, which 
requires the inclusion of pairing interactions, the number of solutions of the full KKDF 
model is twice as great as the number of physical states being described. Hitherto, the major 
technical difficulty (and consumption of cpu time) of this method has been the application 
of a criterion to select the physical solutions. For the ground state problem there is the well- 
known property of the BCS theory that the physical solutions (quasiparticles) correspond 
to positive energies and the unphysical ones to negative energies. In the KKDF model the 
strategy is to ignore initially rotational excitation energies so as to collapse each band to a 
single degenerate state to which the ground state criterion can be applied. We then step up 
the excitation energies, returning them finally to their full values; at each step we select the 
physical solutions by a projection technique described in our cited work, that involves an 
extension of the techniques introduced by Donau and Frauendorf. 

Another way of stating the problem that is directly related to the traditional core-particle 
model is to remark that although only half of the solutions of the KKDF model are related to 
physics, the full set of solutions is necessary for mathematical completeness. The solutions 
of our equations at full excitation can be expanded in terms of the complete set generated 
at zero excitation, but this expansion will involve both physical and unphysical states of the 
latter limit. In the conventional core-particle model it is assumed that the physical states 
of the actual problem are well approximated by a superposition of the physical solutions at 
zero excitation. It follows from this that it suffices to solve a single eigenvalue problem for 
the problem of actual interest rather than having to solve a sequence of such problems. 

The third purpose of this paper is to carry through several illustrative calculations using 
both the KKDF model and the approximation to it just described, in order to assess the 
validity of the latter. In the course of rethinking our algorithms in preparation for this study, 
we have discovered a method of simplifying the full calculation to a sufficient extent that 
much of the advantage of technical simplicity of the core-particle limit has been wiped out. 
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We shall also describe this new development. 

We start in Sees. II and III with a review of the fundamental equations of the Kerman- 
Klein method, in order to introduce some improvements in notation and presentation, as 
well as to correct some phase errors made previously in the formulas for transition matrix 
elements. In Sec. IV we transform our equations (without approximation) to a description 
in terms of an intrinsic frame of reference. Starting from these equations, the definition and 
formulation of the strong coupling core-particle model in its usual form in the intrinsic system 
is given in Sec. V. It is explained in Sec. VI that an equivalent and possibly more effective 
version of this limit can perfectly well be carried out in the laboratory system. Turning 
to applications, our new algorithm is described in Sec. VII and then applied together with 
the standard core-particle model to some illustrative cases in Sec. VIII. Concluding remarks 
are presented in Sec. IX. Two appendices provide some technical details of the derivation 
carried out in Sec. IV. 

II. FUNDAMENTAL EQUATIONS OF THE KERMAN-KLEIN METHOD FOR 

ODD NUCLEI 

We start with a shell-model Hamiltonian of the form 

a 

+jE E F acdh {L)B{ ML {ac)B LML {db) 

abed LM L 

+ E G a UL)A{ ML (ab)A LML (cd). (2.1) 

abed MLm 

Here h a are the spherical single-particle energies referred to the nearest closed shell, a refers 
to the standard set of single-particle quantum numbers, including in particular the pair 
{j a , m a) and a refers to the same set with m a omitted. B\ Mh is the particle-hole multipole 
operator, 

B LM L ( ab ) = E s p (j a m a j b -m b \LM L )a ] a ap 

m a m b 

= {-iy« + i»- ML+1 B L _ ML {ba), (2.2) 
and Al LMl is the particle-particle multipole operator, 

A LM L ( ab )= E 0>V76 - m b \LM L )a ] a a ] p, (2.3) 

m a m b 

where (jimij2m2\jm) is a Clebsch-Gordon (CG) coefficient, s a = (— l)^ - " 1 ", and a bar 
indicates reversal of the sign of the magnetic quantum number. The coefficients F are the 
particle-hole matrix elements, 

Faedb(L) = s y Sfs(j a m a j c - m c \LM L )(j d m d j b - m b \LM L ) V aPl& , (2.4) 

m's 

which satisfies the relation 
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F acdb (L) = (-iy« + i» + ^F bdca (L), 



(2.5) 



and G the particle-particle matrix elements 

Gabcd(L) = ^2(j a m a j b - m b \LM L )(j c m c j d - m d \LM L )V a ^ s , (2-6) 

m's 

which satisfies the conditions 

G acdb (L) = (-l) j " +j c- L+1 G cadb 

= (-iy b+]d - L+1 G acbd . (2.7) 

Our initial task is to obtain equations for the states and energies of an odd nucleus 
assuming that properties of immediately neighboring even nuclei are known. The states of 
the odd nucleus (particle number A) are designated as \Jfj,u), where v denotes all quantum 
numbers besides the angular momentum J and its projection ji. The states of the neighboring 
even nuclei with particle numbers (A±l) are written, in a parallel notation, as \IMn(A±l). 
The corresponding eigenvalues are Ej v and e\^ ±1 \ respectively. The operator equations of 
motion (EOM) are obtained by forming commutators between the single-fermion operators 
and the Hamiltonian, 

[a a , H] = h' a a a 

+ 9 E E Si(j a m a j c - m c \LM)F acdb (L)a 1 B LM {db) 

hd-y LM 

+ \ E E0'«™«Jc - m c \LM)G acbd (L)a\A LM (bd), (2.8) 

1 hdy LM 
[4> H] = ~ h Wa 

~9 E E s i(ia - m a j c m c \LM)F acdb (L)B[ M (db)a\ 

1 hd-y LM 

~\ E E0'« - m a3c m c \LM)G achd {L)A LM {bd)a 1 . (2.9) 

^ hd-y LM 

Here 

K = h a -WF^(L)^±± (2.10) 

4 Lj c Z Ja 1 

are modified single-particle energies. 

The matrix elements of these equations provide expressions that determine the single- 
particle coefficients of fractional parentage, 

Vj^(a; IMn) = (J /iu\a a \I Mn(A + 1)), (2.11) 
Uj^(a; IMn) = (Jfiv\aL\IMn(A - 1)). (2.12) 

To find equations for these quantities, we form the necessary matrix elements of the EOM 
and evaluate the interaction terms by inserting the completeness relation between the single- 
fermion operators and the multipole or pair operators. In order to obtain equations that 
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are expressed completely by means of the amplitudes defined in Eqs. ( |2.11|) and ( |2.12| ), it 
is necessary to interchange the order of the single-fermion operator and the pair operator 
in the interaction terms of ( |2.9| ). This leads to further contributions to the single-particle 
energy in this equation, in that h' a is replaced by h", with 



h" = h' 



V 2L + 1 



Li 



+ l -F 

~ 2 nc<" 



)■ 



(2.13) 



In terms of a convenient and physically meaningful set of energy differences and sets of 
multipole fields and pairing fields defined below, we obtain generalized matrix equations of 
the Hartree-Bogoliubov form 



£ Ju V JpLV {a-IMn) 
£ Jv U Jlxv {a;IMn) 



{e' + u {A+l) + T( A+ V) aIMn , yI , M , n ,Vj^(r, I'M'n') 

+ ^aIMn,-yI'M'n'Uj llu ( / y; I M'n), 

(-e" + ^ A - 1] - r( A -V% IMn ^, M , n ,Uj^(r, I'M'n' 

~ ^alMn^I' M'n'Vj ^vil'i I M Tl). 



Here 



-alMn^I'M'ri 



-E Ju + ^(^ A+1) + E { 

$a~f$irfiMM'5nn'(h' a — \a), 



I 



-E, 



(Ail) 

alNn^I'M'ri 



r (A±l) 

1 aIMn,yI'M'n' 



*aIMn,-yI' M'n' 



Furthermore Eq A±1 ^ refer to the ground state energies of the neighboring even nuclei, the 
matrix elements of are derived from those of ( |2.20| ) simply by the replacement of the 
operator B by B\ and the matrix elements of A' are similarly derived from those of A by 
the replacement of A by together with the interchange A ± 1 — > A =p 1. Finally e a is 
obtained from e a by the replacement of h a by K' a . 

In order to specify a scale for the solutions, we take a suitable matrix element of the 
summed anticommutator, 



(Ail) 



E 



S aJ Sll'5MM'Snn' (Ej n 

- J2 s i(jam a j c - m c \LM L )F acdb (L), 

^ LM L bd 

(I'M'n'(A±l)\B LML (db)\IMn(A±l)) 
I 

- z2(ja m ajc ~ m c \LM L )G acdb (L) 

1 LM L bd 

(I' M'n' (A - l)\A LML {db)\IMn{A + 1)). 



(2.14) 
(2.15) 

(2.16) 
(2.17) 

(2.18) 

(2.19) 

(2.20) 
(2.21) 



n 



a 

E(2ja + 1). 



(2.22) 
(2.23) 



J a 



We thus find 
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O E l\Uj^(a;IMn)\ 2 + \Vj^(a;IMn)\ 2 = l. (2.24) 

All of the above equations are still exact and are not necessarily restricted to deformed 
nuclei. In order to do physics, however, we shall have to impose restrictions on the number 
and nature of the core states included in any application, as well as on the size of the 
single-particle space. 



III. MATRIX ELEMENTS OF SINGLE-PARTICLE TRANSITION OPERATORS 



We next apply the formalism to the computation of matrix elements of single-particle 
tensor operators, T LMl , that we write in the form 

Tlm l = E Vy a /3 a 7- (3-1) 

Pi 

The notation is such that the quantities t a p include a product of matrix elements of single- 
particle operators and of associated coupling strengths (charges, gyromagnetic ratios, etc.) 
We wish to calculate the matrix element {J' h'v'\T LMl \J H>v) ■ To carry through the cal- 
culation, we substitute for the ket a formally exact expression in terms of the action of 
single-particle operators on the states of the core, 

UH = ^ E [Uj,u(a,IMK)al\IMK) 

^ aJMK 



+V Jflu (a,IMK)a a \IMK}}, (3.2) 

where an underline identifies the lighter of the two cores and an overline the heavier one. 
By using the commutation relations and completeness, this leads to the following expression 
for the transition element: 

(J'^'u'\T LMl \J^) = 1 £ [Uj^AaJ'M'K^Uj^aJMK) 

^ a,IMK,I'M'K' 

x( I'M'K '\T t .m l \ IMK ) 



+ [Vj^ lv i(aJ'M'K')Vj^(aJMK){I'M'K'\T LML \IMK) 
+ k E ^[^ V («,M)C/ Jf (7,/MK) 

^ 01,1,1 MK 

-Vj^(a, IMK)Vj llx , v ,{l, IMK)\. (3.3) 

This is now evaluated by use of the Wigner-Eckart theorem with the following definitions 
of the reduced matrix elements: 

(J'pt'u'\T LML \J^) = - fJt\LM L ){J'i/\\T L \\Ju), (3.4) 

(I'M'K'\T LMl \IMK) = \ ' {I'M' I - M\LM L ) 

v 2L + 1 



x(I'K'\\T L \\IK), (3.5) 
V2L + 1 



( — \\3c-m c 

*t»7 = /»,-,-, Uam a j c - m c \LM L )t ac , (3.6) 
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Vj^aJMK) = (^J- (IMJ - fi\ Ja m a )v Ju (aIK), (3.7) 
Uj^(a,IMK) = y (IMJ-^\ Ja m a )u Ju (aIK). (3.8) 

V + 1 

With the help of these definitions, we obtain the formula for the reduced matrix element 
that is utilized in the KKDF model. 



J 1 alKI'K' I 



/ /' L 1 

■/ Ja J 

x [ Mj ,(aJK) Wjv (a/'iT) \T L \\IK) 



+v Jv {aIK)vj, v , (aI'K')(I'K'\\T L \\ IK)} 
+ ^ E U(-l) ia+I+J+i | j j j } uMaIK) Ujl/ (cIK) 

+{ _ 1)ja+I+J+1 1 Jc £ | Vjv{aIK)vj , AcIK) y (3.9) 

This is, with some phase corrections, the formula that was derived in a previous work. 



IV. TRANSFORMATION TO INTRINSIC SYSTEM FOR AXIAL CASE 

We have described previously [l|-|3j several applications of the formalism reviewed in the 
preceding sections to strongly deformed nuclei. Some of the results, together with some 
additional calculations, will be used as the basis for a numerical study of the relation of 
the method of this paper to the traditional strong coupling core-particle model. As will 
be explained in Sec. VI, this relation can be studied using the formalism already at hand 
(theory expressed in the "laboratory" system of coordinates); in fact it turned out to be 
economical for us to carry out all numerical work from this standpoint. Nevertheless, in the 
following two sections we shall undertake to develop the connection between our method and 
the way such calculations are normally presented in the intrinsic system. Our justification 
for this digression is that whenever we have presented a public account of our previous work 
in this field, one question invariably raised was precisely this connection. In what follows, 
we shall answer the question raised in two steps. In the first, we shall derive a form of our 
equations in the intrinsic coordinate system that is fully equivalent to the theory described 
above. Second we shall show that the conventional core-particle approach involves a further 
specialization of this general result and examine this limiting case in some theoretical detail. 

For illustrative purposes, we take a model of the even (core) nuclei that consists of 
the ground-state band \IMK = 0) = \IM) and a finite number of positive parity excited 
bands \IMKn). For the remainder of this section the symbol n will be suppressed. We are 
thus assuming that the eigenstates of the even nuclei have axial symmetry and that their 
eigenstates can be assigned a definite value of K, the component of the angular momentum 
along the figure axis. This assumption is reasonable as long as the states of the same angular 
momentum belonging to different bands are well-separated in energy. 
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We first use rotational invariance to study the structure of the amplitudes V and U 
defined in Eqs. (|2.11| ) and (|2.12|) , respectively. For this purpose we introduce a complete set 



of states \R) localized in the Euler angles, R = (aflj) and write 

\IMK) = J dR\R)(R\IMK) 



27 + 1\* 
8vr 2 



J dR\R)D$ K (R) (4.1) 



The identification of a scalar product of many-body states with the Wigner D function is 



part of the definition of the model. When (4.1) is substituted into the definition of V, and 



use is made of the definitions to follow, we are thereby led to the study of an amplitude such 

as 



(Jfiu\a a \R) = (J^v\U{R)U- l {R)a a U{R)\Q) 

= ^(JfjLulUlJfjL^iJu'ulU^aJJlO) 



E D$ m (R)D^{R)xw v VaK*){-lY-^, (4-2) 



where U(R) is a unitary rotation operator defined by the value of R. The previous ma- 
nipulations have utilized the following relations and definitions (of which the first two are 
standard): 

(JK\U(R)\JM) = D%>(R), (4.3) 

U-\R)a 3m U(R) = J2 a jK D<£>{R), (4.4) 

K 

(JHSJO) = (-iy +m XJ,u(jm) (4.5) 



The introduction of the phase in ( |4.5|) simplifies the structure of the transformed equations 
of motion given below. 

With the help of the integral of a product of three D functions and the application of 
standard symmetry properties of CG coefficients, we find 



57T 2 



V Jfa ,(a; IMK) = £ J ——-(-1) J ^(IMJ - /i\j a m a ) 

K a V ^Ja + J- 

x(JK - K a3a K a \IK){-iy^ X JK-nAja^a). (4.6) 

A similar analysis carried out for the amplitude U yields the result 



57T 2 



Uj^ia; IMK) = £ ^__(-l)^+^-««+^+^(/MJ - fi\ Ja m a ) 



X(JK - K a jaKa\IK)^jK-Kau(jaKa), (4.7) 
0J^(ja« a ) = (^H a },-«al°)- ( 4 - 8 ) 

Starting from Eqs. ( |2.14|) and ( |2.15| ) and utilizing the forms ( |4.6|) and ( f4.7|) , we next 
derive equations satisfied by the amplitudes x an d 4>- The technique is to eliminate the CG 
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coefficients that occur in (|4.6| ) and (|4.7|) by multiplying by {I M J — n\j a m a ){J K — K a j a n a \I K) 
and by the reciprocal of the factors pre-multiplying these CG coefficients in the one or the 
other of these equations, summing over M, \i and /, and using standard formulas of angular 
momentum algebra. Some details are provided in Appendix A. In the equations to follow, 
the quantities that appear for the first time are defined by the equations 



K(m, K\j, J) = \J(j + m)(j - m + 1) 



xJ(J - K + m)(J + K -m+1), (4.9) 



[I'M'K'\B[ M (db)\IMK) = q < £$(db) ] 



21+1 



x (IMLM L \I'M')(IKLK' - K\I'K'), (4.10) 

(I'M'K'\A{ ML (db)\IMK) = ASVfe)^ 

x (IMLM L \I'M')(IKLK' - K\I'K'), (4.11) 

= E%**> + -^[1(1 + 1) - n (4-12) 
21 K 

Of these equations, the quantity 1Z is recognized as arising from the matrix elements of 
the Coriolis coupling and the remaining equations are expressions valid for the axial rotor 
model for matrix elements of transition operators (see further below) and excitation energies. 
These expressions constitute definitions of the intrinsic multipole moments q, of the intrinsic 
pairing moments A, of the band-head energies Ek, and of the moments of inertia Ik- 
The resulting equations (with partial suppression of the index v) are 

£j»Xj,K- Ka (jaKa) = K + E { t l) + —^njlAJ +l)~K 2 

+ja(ja + 1) + 2K a (K - K a )]}xj,K- Ka (jaK a ) 
+ (A+l) 7 ^"' K \j*' J)XJ,K-K a +l{j a K a ~ 1) 

2Z K 

+ {A +l) n (~ Ka ^ ~ K \j^ J)XJ,K-Ka-l{jaKa + 1) 

2Z K 

bcdn c K' L 

X{jc ~ K c j a K a \LK - K')xj,K-K a {jcK c ) 

+ E ^(-l) Jc+Ka+L G QC ,6(^)AK ) (d6) 

bcd,K c K' L 

x (jc ~ K c j a n a \LK - K')(j)j jK _ Ka (j c K c ), (4.13) 

£jAj,K- Ka (jaKa) = {"< + + [J {J + 1) ~ & 

21 K 

+ja(ja + 1) + 2K a (K - K a )}}4>j,K-K a {jaKa) 
+ (A-l) ^(^' K \3a, J)(pJ,K-K a +l(jaHa ~ 1) 

2X K 
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+ ~ K \3a, J)<Pj,K-K a -l{j a K a + 1) 

2X K 

- E ^(-i^^W^C*) 

bcd,K c K' L 

X(jc - KcjaKa\LK - K')<fij iK - Ka (j c K c ) 

+ E ^ rf (i)AS(*)(-)^- 

X (j c - K c j a K a \LK - K')xj,K~K a {jc^c)- (4.14) 

In these expressions, we have deliberately chosen, for conciseness of expression, not to do 
the sum on k c , where the value k c = K' — K — K a is imposed by the resident CG coefficient. 
Relations ( |4.10| ) and ( |4.11 ), which have been used in all our previous applications, are 



approximate, and therefore require further discussion. For example, Eq. (|4.10| ) follows as 
the value of the first term of the operator expression 

Bl ML W= E Y< P i>M>K>q i X ) (db) 

Xi,...,X p 

i' Xi ...r x jp IMK . (4.i5) 

Here the P are the projection operators for the specified band members, I' x is a spherical 
tensor component of the intrinsic angular momentum, and the braces imply a symmetrized 
expression. Assuming that the connected bands have the same parity, p is even for even 
electric multipoles and odd magnetic multipoles and odd for odd electric multipoles and even 
magnetic multipoles. If the connected bands have opposite parity, there is a corresponding 
relation. The form of (|4.15 ) is a consequence of the assumption that B must be a tensor 



operator of appropriate rank in the Hilbert space of the axial rotor. The further assumption 
that we can limit ourselves to the first term is that for the states of interest the rotor is almost 
rigid, as is true for the low-lying states of strongly deformed nuclei. The corresponding 
expression for the pairing operator requires only the replacements 

f^-Afeg, (4.16) 

and the realization that the projection operators to the left and to the right refer to different 
cores. 

The inclusion of odd multipole or pairing interactions requires that, minimally, we choose 
p — 1. The evaluation of such a multipole term is carried out in Appendix B. 



V. CORE-PARTICLE COUPLING MODEL 



A. Spectra 

For further development, we specialize the formulas of the previous section to the con- 
ventional monopole pairing plus quadrupole-quadrupole model and confine our attention 
initially to the special case that we include only the ground-state band of the neighboring 
even cores. (The general case will be considered subsequently.) We also assume that we are 
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treating well-deformed nuclei and ignore number conservation. For L = pairing we have 
in the limit of a constant pairing matrix element 



]T G aabb A§; 0) (bb) = 2A a J2j a + 1 



= 2A v /2j a + 1. (5.1) 

For the quadrupole interaction, we write 

F a bcd(2) = —^FabFdc, (5.2) 

^2F db qg'°\db) = Q . (5.3) 
bd 

Because we are dealing with a K = band, axial symmetry implies that n a = k c = k, and 
the quadrupole potential becomes 

K c = -\K 2 F ac Q (-iy^( 3c - Kj aK \20). (5.4) 

The potential V is symmetric provided we choose 

F ca = (-l) j ^ +1 F ac , (5.5) 

which is consistent with ( |2.5|) . 

We next study the limit of our equations found by introducing the simplifications made 
above and also neglecting the core excitation energies. The resulting equations do not depend 
on the total angular momentum, and we thus set (with n a = k) 

XJ-Ku'c^) * "X.KC1 

0J,-«O*c«) ~» 0/CC, 

-> £ K r- (5.6) 

Evidently k is the component of the quasi-particle angular momentum along the axis of sym- 
metry, and r resolves degeneracies in the values of k. In the limit considered our equations 
thus reduce to a Hartree-Bogoliubov set 

SktXkcl = taXna + ^IcXna ~ A0 Ka , (5.7) 
S KT (f) K a = -e a (j) Ka - VacCpKa ~ A^Ka- (5.8) 

^From now on we set e' a = e" = e a . 

These equations are solved by introducing the unitary transformation that diagonalizes 
the single-particle Hamiltonian 

K c = ej ac + V* , (5.9) 

namely, 
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T 

v kt = J2 AK c?Xkc (5.10) 

c 

^ A ar A bT = $ab, 

T 

a 

£A^ a V^ = ^rr'. (5.11) 



We thus obtain a standard set of BCS equations 

£ktV KT = €-KT V KT ^KT) (5-12) 

£-k,tU kt = Crt^KT ^KT) (5.13) 

with the usual solutions 

S KT = ±v/e2 T + A2, (5.14) 
where corresponding to the plus sign, we have the physical solutions 

*-={£)• (5 ' 15) 

and to the minus sign the unphysical solutions 

We have reviewed this familiar material because of its importance in the definition of the 
standard core-particle model. 

We have now laid the groundwork for the solution of the full equations of motion Q4.13 ) 



and (|4.14 ). For this general solution the notational change contained in ( p.6[ ) is generalized 
to 

XJ, — K.{jc^) * XjKC) 
0J,- K (j» -> (f>jKC, 

£jv — > £jkt- (5-17) 
Introducing again the transformation that diagonalizes the single-particle Hamiltonian Ti, K , 

T 

XjKT = Y, A arXJ K a, (5.18) 

a 

with a corresponding transformation for (ft, the equations of motion become 
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^JktXJkt &kt~X,J kt A®Jkt 

+ J2Ui, K , T , X J K >T>, (5.19) 

k't' 

£jkt4>Jkt = —£kt4>Jkt — A-XJkt 

+ E^v^«v, (5.20) 

k't' 

and the non-vanishing matrix elements of U that occur in these equations (that reinstate 
the angular momentum and include the Coriolis coupling) are 

2TU J „, = J2 A* ar [J(J + 1) + j a (j a + 1) - 2K 2 ]A ar ,, 

a 

2iui tK _ lT , = KM*, o|j a , J)A aTl , 

a 

2XUi tK+lT , = KM~^ 0|j a , J)A aT , (5.21) 

a 

We recall that the quantities 1Z are defined in ([4.9|) . 

At this point we introduce the defining approximation for the core-particle coupling 
model, first defining 

< 5 - 22 > 

and then assuming that \l/ can be expanded in terms of the physical solutions of Eqs. ( |5.7|) 
and (|5.8|), namely, 

Vjkt = Cj KT iP KT . (5.23) 

By contrast, the exact expression must be the sum of a physical and an unphysical solution. 
If we include the latter, we have, in fact, returned to the KKDF model and to its basic 
technical problem of selecting physical solutions. Although, as we shall see later, we have 
found a simplified method to handle this problem, compared to the approach used in earlier 
work, it remains of interest to know when the traditional core-particle model is valid. 

With the help of ( |5.23|) , Eqs. ( |5.19|) and (|5.20|) can be reduced to a form of the core- 



particle coupling equations ready for final numerical study, namely, 

£jktCj kt = £ KT Cj KT + ^2 Wk TK / t >C "j k > t > (5-24) 

k't' 

W kt ,k't' = Tp KT U^ TiK , T ,ip K > T/ . (5.25) 

This is a standard diagonalization problem with the "correct" number of solutions. These 
solutions will be compared with the exact solution of the corresponding KKDF equations. 

We consider next the general case defined in the theoretical formulation of the previous 
section, with multiple bands in the core nuclei, but with the maintenance of axial symmetry. 
Though not really necessary, it makes sound physical sense to proceed as follows: We lean 
on the fact that the interband quadrupole transitions are weak compared to intraband 
transitions. Thus we shall first ignore the terms associated with these transitions as well 
as the perturbation associated with finite excitation energy above the band-head. What 



13 



remains is a Hartree-Bogoliubov approximation for excited bands. Next we add the "Coriolis 
coupling" and thus obtain a series of bands in close analogy with our treatment of Coriolis 
coupling for the ground-state band. Finally, we introduce the coupling arising from interband 
transitions in the cores. 

In fact, it is hardly necessary to give many details of the previous steps. All we need 
is an enhanced notation. Instead of the ground state band, we consider a band Ka, where 
00 is the ground state band, 01 the beta band, 20 the gamma band, etc. Now to all the 
quantities defined above, such as A, Q , V„ c , %kci 4>kci etc., we add a superscript {Ka). Thus 
after transformation by the matrix A^ Ka the excited state HB equations become 

- E K °)v« T ° = e K K °v K J - A K °u«°, (5.26) 
(£«° - E K °)u K K : = -egu** - A K °v«°, (5.27) 

with the solutions 



where the first term on the right hand side is clearly the band-head energy. The remainder 
of the calculation also parallels that made for the case of the ground-state band. The only 
quantities requiring more than a notational change are the matrix elements of the operator 
U defined in Eq. (|5.21|) . The necessary emendations can be read off directly form the core- 
particle equations ( |4.13| ) and ( |4.14| ). 

Thus we have specified a procedure for deriving a set of state vectors HfjJ and associated 
energies £f°, where we have amalgamated the pair of quantum numbers kt into the symbol 
v. We have taken account of all terms in the effective Hamiltonian except for the interband 
multipole fields. To finally include the latter, we write 

V = r 3 (V d + V od ), (5.29) 

where t% is the usual Pauli matrix, d refers to the intraband parts of the multipole field and 
od to the interband parts. It remains to take into account only the latter piece. This is done 
by a final mixing 

®Jp = E A J wKa*%, (5-30) 

vKcr 

where the mixing coefficients are determined by the conditions 

c aJ _ fKaU 



+ X/ J~l K(TV , K , <J ,A J p y K i a i, (5.31) 

v'K'o' 

Hk^k>«> = ^r 3 V od ^y. (5.32) 



In the last two sections, we have derived the conventional form of the core-particle 
coupling theory from the KKDF formalism. In fact the equations derived in the first of these 
sections were exact, i. e., completely equivalent to those of KKDF, indeed only their form in 
the "intrinsic frame". The core-particle coupling model as customarily presented involves, 
as described, an additional approximation in the solution of these equations. Indeed, the 
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essence of the model lies in this approximation rather than in whether calculations are carried 
out in the intrinsic system as described above or in the laboratory system as is done in the 
full application of the KKDF method. In the next full section, we shall record the form of 
the core-particle approximation in the intrinsic system, coupling theory in the laboratory 
frame. 



B. Core-particle coupling model: transitions 

Here we shall only indicate the step involving the exact transformation of ( |3.9|) into 
an expression referring to the intrinsic system. Since we shall not utilize this version of 
the formalism, we leave the further transformation by the introduction of the approximate 
solutions developed in the preceding subsection as an exercise for the reader. This step is to 
introduce values for the reduced matrix elements on the right hand side of (3.9) and to carry 
out the summations over / and I' in order to reach a formula appropriate to the core-particle 
coupling model. By comparing (|3.4j) with (|4.10|) , ( p.7|) with ( fOD , and ( |3.8|) with ( [4.7|) , we 
can read off the formulas 

(I'K'\\T L \\IK) = V2TTTq { K$(IKLK' - K\I'K'), (5.33) 
v Jv (alK) = y/^{-l) ja+Ka {JK - K a3a K a \IK) 

K a 

*XJK- Ka (jaKa), (5.34) 

u Jv {aIK) = £ V8^(-iy*- Ka (JK - Kj aKa \IK) 

ft a 

X(f)jK-K a UaKa)- (5.35) 

Carrying out the summations over / and we are led to the equation 

(J'v'\\T L \\Jv) = ^ Y ; 1 (JK — K a LK' — K\J'K' — Kg) 



X [(pJK~K a {jaKa)4>J>K>-K a (j 



+ ~FT E ^c[(-l) ic+Kc (ja - K a j c K c \Ln c - K a ){JK - K C LK C - K a \J'K - K a 



acK 



/ (2L + 1)(2J' + 1) 



+ (-l) ic+Kc+J+J ' +L 0'a - K a j c K c \LK c - K a ){J'K - K C LK C - K a \JK - K a ) 

X 7= w XJK- Ka (jaK a )XJ>K~K c (jc^c) ■ (5.36) 

/ (2L + 1)(2J + 1) 



VI. CORE-PARTICLE COUPLING MODEL IN LABORATORY FRAME 

We show here that the core-particle coupling model can be formulated just as conve- 
niently in the laboratory frame of reference as in the intrinsic frame. We start with the 
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fundamental matrix equations of motion, Eqs. (|2.14j) and fl2.15|) and reduce them by ap- 
plication of the Wigner-Eckart theorem. By means of Eqs. (|3.7|) and (|3.8|) , Eqs. ([4.10 ) 
and ( |4.11| ), and standard angular momentum algebra, we find the equations (assuming that 
K' — K and L are even, as is the case for the specific model considered in the body of this 
paper), 

E Jv v Jv {aIK) = (e a + u IK )v Jv (aIK) 

+ Y, T(aIK,cI'K')v Ju (cI'K') 

el'K' 

+ ]T A(aIK,cI'K')u Jv (cI'K'), (6.1) 

el'K' 

£j u u Ju {aIK) = {-e a + u) IK )ujv{aIK) 

- T{aIK,cI'K')u Jv (cI'K') 

el'K' 

+ Y^IK,cI'K')vj v (cI'K'), (6.2) 

el'K' 

T(aIK,cI'K') = WF a *u>(L)q%$(db)yf(2L+l)(2I+l) 

Z Lbd 

x (-lyo+'+J | 3* J j: L j | (IKLK' - K\I'K'), (6.3) 
A{aIK,cI'K') = W G acdb (L)A^(db) y /(2L + l){21 + 1) 

^ Lbd 

x (-lyo+t+J I 3* j j: L \ (IKLK 1 - K\I'K'). (6.4) 



In Eqs. ( |6.1| ) and ( |6.2D we have set e' a 

We introduce a condensed notation for (|6.1|) and (|6.2| ), writing them in the operator form 

e Jv ^ Ju = fc^ Jv + uVj u , (6.5) 

*=(l r _ e A _ r ), (6-6) 

We solve these equations in the approximation that forms part of the definition of the 
core-particle model. Again we consider first the simplest case where only the ground state 
band of the cores is included. The extension to excited bands and interband coupling can be 
dealt with in analogy to the treatment described for the intrinsic system. Let the physical 
solutions of (|6.5|) with u = be designated as ^jj, with corresponding energies 8f}. Here 
the symbol v abbreviates the set (act). We approximate the solutions of the full equation 
by the expansion 

* J* = (6.8) 



The introduction of this expansion into ( |6.5|) leads immediately to the standard eigenvalue 
problem 
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£jvCt u i — Zjv'Ctu' + v iiC j vv ii , (6-9) 

u" 

U J V , V ,, =¥$0,*®,. (6.10) 

This equation is to be compared with Eq. ( |5.24j ), to which it is equivalent as long as u has 
the form assumed in the derivation of the latter. In fact, ( |6.9|) has an advantage in the 
case that the excitation spectrum is not conveniently expressed in algebraic form, but its 
numerical values are known from experiment. 

We can extend the theory to include multiple bands in the core nuclei. We use the 
labels Ka to distinguish the different bands and now take as a zeroth approximation the 
coupling of the odd particle to a single one of these bands. The theory is, to start with, the 
same as that described above except that we must distinguish the results for the various 
cores, and this is done by a superscript Ka. In so far as the multipole fields and pairing 
fields for the band Ka are almost equal to those for the ground band, the energies Sfj 
are almost independent of Ka. We prefer to lift this degeneracy by shifting each of these 
energies by E Ka , the band head energy, and redefining uj to be the excitation energy above 
the band head in each case. The step that follows is to introduce the mixing due to the core 
excitations and again only the change in notation already specified is necessary to record 
the equations that generalize Q6.9Q and fl6.10|) . 



The final step is to include the further mixing due to interband multipole fields (assuming 
that such mixing for the pairing fields can be neglected). For this purpose, we decompose V 
into an intraband piece (subscript d) and an interband part (subscript od), the latter having 
so far been neglected, according to the equation 

f = r 3 (f d + f od ), (6.11) 

where r 3 is the usual Pauli matrix. The perturbation previously neglected is dealt with by 
the expansion 

0;p = E ^pW^> ( 6 - 12 ) 

vKo 

where the mixing coefficients are determined by the conditions 

C t>J _ cKcT'p i J 
°Jp u p,vKo — °Jv u p,vKa 



+ E GvKcr,v , K'a''DpyK'a'i (6.13) 
v'K'a' 

QiwK'a' = *%°T 3 t od *%?. (6.14) 



VII. IMPROVED ALGORITHM 

The main source of difficulty perceived in the solution of the KKDF equations is that 
the set of solutions is over-complete by a factor of two. This is a consequence of the fact 
that the basis states form an over-complete (and, consequently, non-orthogonal set). Thus 
half of the states found by solving the EOM are not physical and have to be identified and 
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removed. The technique previously used to perform this task has now been understood to 
be unnecessarily complicated. 

In the previous approach [|1],[7| the Hamiltonian is first decomposed into symmetric and 
anti-symmetric parts with respect to particle- hole conjugation. If only the anti-symmetric 
part is diagonalized, then for every positive energy eigenvalue there is a negative partner. 
From the BCS theory we know that the positive eigenvalues are the physical solutions and the 
negative eigenvalues the non-physical ones. Then the symmetric part is turned on "slowly" 
and at every step the physical solutions are identified using a projection operator built from 
the wavefunctions of the previous step. Since the equations of motion have to be solved at 
each step, the time needed to perform the calculation is correspondingly longer than for a 
single diagonalization. (In most applications a typical number of steps is 5.) 

A simpler and quicker approach has now been identified. Since the problem decomposes 
into subproblems involving states of a fixed angular momentum, we can invoke the no- 
crossing theorem. This means that the relative order in energy of the physical and of the 
non-physical states does not change as we turn on the symmetric part of the Hamiltonian. 
If the lower half of the states (negative in particular) are the unphysical ones in the BCS 
limit, then at the physical limit where the full Hamiltonian is used, the lower half of the 
states are again the unphysical ones. Consequently, we need only to solve the equations 
of motion at the two limits, the BCS limit and the full Hamiltonian limit. These remarks 
about the technique of solution apply not only to the strong coupling examples studied in 
the next section, but also to less straightforward applications of the KKDF method. 



VIII. APPLICATIONS 

We illustrate the remarks of the previous sections with applications to a pair of well- 
deformed nuclei. The first application is to the nucleus 157 Gd, which we have studied pre- 
viously 0,3. 157 Gd is a well deformed nucleus and thus suitable for application of the 
strong coupling core-particle model. To recall a few details, we used a large single-particle 
space (including all states from 5 major shells). The energies and matrix elements of these 
single-particle levels were calculated using the Woods-Saxon potential. The odd neutron is 
coupled to the cores 156 Gd and 158 Gd, which are represented not only by their ground bands, 
but also by several excited bands, as was found necessary to fit all the observed bands of 
157 Gd. The core excitation energies, cuj, were given by phenomenological formulas tuned 
to experiment. In the same way as in the previous papers, the strength of the quadrupole 
field is treated as a free parameter and the values of the single-particle energies found from 
Woods-Saxon calculations are allowed to vary by ±5%. First we solved the EOM problem of 
the full KKDF model and fixed the strength of the quadrupole force and the single-particle 
energies in order to achieve the best fit. Then we solved the EOM for the core-particle model 
as described in Sec. |VT], using the same parameters. The results are show in Fig. ([[J). We 
can see from the figure that the two models give very similar results. In Fig. @ we show 
the result of the B(E2) calculations. Again it is clear that the two models give very similar 
results. 

The second application was to the proton spectrum of 157 Tb, with 156 Gd and 158 Dy cores. 
We used the same method as described above and the results are shown in Fig. (|2|). The 
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conclusion is the same as in the previous application, namely that the two methods give 
very similar results. Observed B(E2) values are too few to allow a meaningful comparison. 

To the extent that the examples chosen are typical, it is apparent that for well-deformed 
nuclei the strong coupling core-particle model gives almost as good results as the full KKDF 
model. We emphasize, however, the greater range of validity of the KKDF model, in partic- 
ular to cases such as transitional nuclei PUl8||, where none of the usual traditional versions 
of the core-particle model is applicable. 



IX. DISCUSSION AND CONCLUDING REMARKS 



In this paper, we have studied a semi-microscopic core-particle coupling theory, the 
KKDF theory, and particularly its relationship to the traditional strong coupling core- 
particle model. The KKDF theory is formulated in the laboratory system of coordinates, 
and as such, can be applied both to the spherical vibrational (weak coupling) and deformed 
rotational (strong coupling) regimes, as well as to transitional cases. A significant portion 
of this paper has been devoted to transforming the KKDF equations from the laboratory to 
the intrinsic system of coordinates, the latter defined only for the well-deformed regime. We 
have pointed out the additional approximation necessary to reduce the KKDF equations to 
those of the usual core-particle limit. We have then applied both the full and the limiting 
model to a few illustrative nuclei and found only small differences in the numerical results. 
This justification is, for our purposes, less significant than it would have been in the past, 
since we have also formulated an improved algorithm that renders the KKDF equations 
essentially as simple to deal with as the defined approximation. 

The reason for the good agreement between the approximate and the complete theory 
obviously expresses the fact that there is little mixing between physical and unphysical states 
as we "turn on" the coupling that is initially suppressed in our approach. This means that 
they stay well separated in energy. We can expect this situation to change for applications 
where there are multiple avoided crossings. 



APPENDIX A: SOME DETAILS OF THE DERIVATION OF THE 
CORE-PARTICLE COUPLING MODEL 



We provide some details of the derivations of Eqs. (|4.13|) and ( |4.14|) . The first terms that 
require special attention are those involving the excitation energy in the even nuclei. We 
immediately do the sum over M, {i. Now consider Eq. ( 4.13 ), where we encounter the term 



(u IK - E K )(Jfi'j a K a \IK) = Ua - ««J/x'|r^[(J + j) 2 - K 2 ]\IK), 

(suppressing mass number). We can replace the combination (J + j) 2 by 
J(J + 1) + j a (j a + 1) + 2(K - K a )n a + j-J+ + j+ J_. 



(A.l) 



(A.2) 



Applying the standard algebra of the raising and lowering operators and shifting the variables 
n a as required for these terms, we thus obtain additional contributions of single-particle type 
as well as the Coriolis coupling. 



19 



We consider next the contributions of the multipole and pairing fields, a calculation that 
requires most of the modest labor involved in the derivation of Eqs. ( 4.131 ) and ( 4.14 ). As 
an example of what is involved, we compute the contribution of the even multipoles to the 
right hand side of (|4.14| ), which we label T(T* : JKj a K a ). Utilizing Eqs. ( ET7D and ( |4.10|) , 
we must evaluate the expression 



T(rt : JKj a K a ) = ^J2(-iy c+ma+Kc - Ka 



(2j + i)(2/ + i; 



\J (2j c + l)(2/' + 1) 
xF acdb (L)q^(db) 

x(j a ~ m a j c m c \LM L )(IMLM L \I , M')(I , M'Jm c - M'\ Jc m c ) 
x(IMJm a - M\ Ja m a ){JK \IK) 

x{IKLK' K \I'K'){JK' - K c j c K c \l'K')<t>j, K ,- Kc {jcKc). (A.3) 

In this equation the sum is over all angular momentum variables not indicated explicitly on 
the left hand side except for m a , which disappears from the final result. 
To evaluate this expression, we first study the partial sum 



5 = EC- 1 



\jc+m a 



I2 ja + 1 



V 2 j c + 1 

x(j a - m a j c m c \Lm c - m a )(IMLm c - m a | I'M') [I'M' Jm c - M'\j c m c ) 
(-l) J '» +/ ' +1 V(2L + l)(2/'+l) 

w J ja jc L 



I' I J 



(IMJm a - M\j a m a ) 



(A.4) 



which can be derived from Edmonds (6.2.7). The sum over M then removes two more CG 
coefficients from ( |A.3| ) The next step is to apply Edmonds (6.2.6) to evaluate the sum over 
leading to a final trivial sum over /. We thus find 



S' = Y.(- i y a+V+J ^ 2L + !)( 2/ + l )i J K - K a3a K a \IK) 



ii 1 



x(IKLK' - K\I'K'){JK' - k c3c k c \I'K') j 3 f ^ j 
(-l)"°-"o+*°+"°(j c - K c j a K a \LK' - K). 



(A.5) 



Equations (|A.4| ) and (|A.5|) are the essential results for the evaluation of ( |A.3|) leading to 
the appropriate term in ( [4.141) . The pairing term in the same equation and the multipole and 
pairing terms in ( |4. 13| ) can be shown (after straightforward transformations for the latter) 
to involve the same basic sums, up to phase factors. 



APPENDIX B: CONTRIBUTION OF ODD MULTIPOLE OPERATORS 



In the main text, we have included in the general core-particle equations ( |4.13| ) and 



(4.14) only contributions from even electric multipole-multipole interactions, assuming that 
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the included bands are all of the same parity. If we wish to include odd electric multipole- 
multipole forces to lowest order, we must replace the matrix element Eq. ( 4.10 ) by the value 
of the matrix element 

(I'M'K'\B{ ML (db)\IMK) 

= q^(db)Y,(l'M'K>\{Dl lLK ,_ K+x J' x }\IMK). (B.l) 

A 

Here I' x , X = ±1,3 are the spherical tensor components of the angular momentum in the 
"intrinsic" system, expressed in terms of the Cartesian components (where they differ) by 
the equations 



J'-i = ^« - »k) = (B.2) 

With the help of the well-known matrix elements of the intrinsic components, Eq. (|B.1| ) 
takes the value 

(I'M'K'\B[ ML {db)\IMK) = q%$(db)(IMLM L \l'M') 
1 



x{ — y=[{IK - ILK' -K + l\I'K')^J(I + K)(I-K + 1) 



+J(P - K')(I> + K'+ \){IKLK' -K + 1\I'K' + 1)] 



+-=[(IK + ILK' -K- 1\I'K')^{I - K)(I + K + 1) 
V & 



\I' + K')(r -K' + 1){IKLK' - K - 1\I'K' - 1)] 
+(K + K')(IKLK' - K\I'K')}. (B.3) 

With these values, we are now in a position to evaluate the contributions of an odd electric 
multipole force to our core-particle coupling equations. We first consider the contributions 
to ( |4.14D . As an example, consider the first term of ( |B.3| ). The calculation parallels that 
described in the previous appendix. The sum ( A.4 ) repeats itself in every case. The sum 
( |A.5| ) is replaced, in general, by different expressions. In the case of the first term of ( |B.3[ ), 
the sum S' is replaced by the sum 



S' X = J2(- 1 ) ja+I ' +J i( 2L + !)( 2/ + l )( JK ~ Ka3aK a \IK) 



x(IK - ILK' -K+ 1\I'K')(JK' - K c j c K c \I'K' 

w I 3a jc L 



V / j^(I + K)(I-K+l). (B.4) 

For S[, the sum over I' can be carried out as before, but after this has been done, instead 
of a final normalization condition for CG coefficients, we encounter the sum 



J2(JK - K a j a K a \I K)(JK - Ka - ljaK a \IK - 1) J (I + K) (I - K + 1 



(J- K + Ka + l)(J + K-k), (B.5) 
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which involves the same "trick" as used in the evaluation of the Coriolis coupling. Of the 
terms arising from ( ]B.3|) , the first and third require the procedure just described, the second 
and fourth a similar procedure in which we interchange the order of the sums on I and 
and the fifth the same calculation as in the previous appendix. We also find that the first 
two terms are equal, as are the third and fourth. 

Altogether, we find for the contribution to Eq. (|4.14j ), the expression 

E (-iy c+Kc F acdb (L)q ( K L $(db) 

bcd,K c K' L 

X [-^=\J (J ~ K + K a + 1)(J+ K - K a )(jc ~ KcjaK a \LK — K' — l)(f)j jK ^ Ka (j c K c ) 
~~^2^ J + K ~ Ka + ~ K + Ka ^ C ~ Kc i aKa \ LK - K ' + l)<f>J,K+l-K a (jcKc) 

~(K + K')(j c - K c]a n a \LK - K')(/>j >K _ Ka (j c K c ). (B.6) 

For conciseness of expression, we have not done the sum over k c . In this form it can be 
shown that the corresponding contribution to Eq. (|4. 13| ) differs only by overall sign and by 
the replacement 
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FIG. 1. Negative parity energy levels for 157 Gd. The circles correspond to the experimental 
values, the solid line to the KKDF model and the dotted line to the core-particle coupling model. 
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FIG. 2. Positive parity energy levels for 157 Tb. The circles correspond to the experimental 
values, the solid line to the KKDF model and the dotted line to the core-particle coupling model. 
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FIG. 3. B(E2) transitions for 157 Gd. Comparison of the KKDF model and the core-particle 
coupling model. The points with error bars are the experimental data, the dashed lines result from 
the core-particle model and the solid lines from the KKDF model 
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